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We examine cluster states transformed by stochastic local operations and classical communication, 
as a resource for deterministic universal computation driven strictly by projective measurements. 
We identity circumstances under which such states in one dimension constitute resources for random- 
length single-qubit rotations, in one case quasi-deterministically (N — U — N states) and in another 
probabilistically (B — U — B states). In contrast to the cluster states, the N — U — N states exhibit 
spin correlation lunctions that decay exponentially with distance, while the B — U — B states can be 
arbitrarily locally pure. A two-dimensional square N — U — N lattice is a universal resource for quasi- 
deterministic measurement-based quantum computation. Measurements on cubic B — U — B states 
yield two-dimensional cluster states with bond defects, whose connectivity exceeds the percolation 
threshold for a critical value of the local purity. 



I. INTRODUCTION 

In the Measurement-Based Quantum Computation 
(MBQC) model [TJ [2], one starts with a highly entan- 
gled many-qubit quantum state called a resource state, 
and processes logical information via single-qubit mea- 
surements on the physical qubits of the resource state. 
In order to compensate for the randomness of the mea- 
surement outcomes, the bases in which measurements are 
performed must be conditioned on the outcomes of previ- 
ous measurements. Proceeding in this way, one can tele- 
port logical quantum information situated on one part 
of the state to another part, but having been subjected 
to some desired unitary transformation. If the basic uni- 
tary transformations that can be applied via single-qubit 
measurements on the resource state generate a set that is 
dense in SU(2), then the resource is said to be universal 
(in the terminology of Ref. [3J , this is the notion of CQ- 
universality, and such a resource state would be called a 
universal state preparator). 

The archetypal family of resource states known to 
be universal for efficient MBQC is the so-called cluster- 
state [HE]- This state is special in several ways: all spin 
correlation functions are strictly nearest-neighbor [H [5] , 
the localizable entanglement between any pair of qubits 
is maximal 01 E], it is the only state (up to local uni- 
taries) on small system sizes that saturates the Tsallis 
and Renyi entropies of entanglement [5] , it cannot be the 
non-degenerate ground state of a two-body spin Hamil- 
tonian [5] , and so on. One might expect that one or 
more of these properties would be necessarily satisfied 
by any universal resource state. This has turned out 
not to be the case; several authors in recent years have 
identified resources that differ materially from the cluster 
states [H El EMU]- 

The newly discovered richness in the landscape of re- 
sources notwithstanding, the property of universality is 
exceedingly rare; not only must a family of universal re- 
source states saturate various measures of entanglement 
in the thermodynamic limit [31 [IS], but the entangle- 



ment with respect to other measures must not be too 
high [16l [17] . Therefore, it is highly unlikely that a ran- 
dom pure state will be universal, so a search for new 
resources must be heavily constrained in order to have a 
reasonable chance of success. 

Recently, a number of new resources [U 02 [101 [18] have 
been proven to be universal by means of reduction to a 
known resource state [TS]. The reduction strategies of 
interest are those composed purely of local operations, 
possibly augmented by classical communication. They 
are typically stochastic, in the sense that the known re- 
source state is smaller than the original state. In other 
words, these resources all appear to be within the equiv- 
alence class of the cluster states under Stochastic Lo- 
cal Operations and Classical Communication (SLOCC). 
The SLOCC-equi valence class of an n-qubit pure state 
is known to be its orbit under GL (2, C)®", the group of 
ri-fold tensor products of two-by-two invertible matrices 
over the complex numbers. In other words, two n-qubit 
states \tp) and \(j>) are equivalent under SLOCC if and 
only if 

\xl>) = S {1] ®S 2 ®---®S n \$) (1) 

where the {Si} € GL(2, C) are invertible, two-by-two 
complex matrices. 

A natural question thus arises: are all universal MBQC 
resources SLOCC-equivalent to the cluster states? More 
precisely, is any n-qubit element of a family of universal 
resource states SLOCC-equivalent to an n-qubit cluster 
state? 

In this paper, we tackle a related question, namely: 
what states in the SLOCC-equivalence class of the two- 
dimensional cluster states are universal for MBQC? It is 
clear, by construction, that each state in this class can 
be stochastically reduced to a cluster state, but what is 
not clear is whether it is possible to compute directly on 
the image of a cluster state under some invertible, lo- 
cal map. We show that there is a restricted subclass of 
invertible local transformations, strictly including the lo- 
cal unitaries, whose image is a set of quasi-deterministic 
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resources for MBQC, where in general the computation 
is of random length and 'repeat-until-success' strategies 
must be employed (c.f. Refs. [5[[T5]). In particular, we 
identify two types of SLOCC operators whose action can 
in certain cases preserve the usefulness of the cluster state 
as a resource. The first type, which we call N-type opera- 
tors, comprises those operators that preserve the relative 
norms of the computational basis states. The second 
type, called B-type operators, are those that preserve 
their orthogonality (i.e. are in a sense basis-preserving). 

In particular, we show that when N-type operators act 
on alternating qubits in a ID cluster state, the state re- 
mains a quasi-deterministic resource for single-qubit ro- 
tations. We refer to such ID states as N— U— N chains. In 
contrast to the cluster state, the number of measurements 
required to implement an arbitrary single-qubit rotation 
with an N — U — N chain is random rather than fixed. 
Furthermore, the state exhibits non-zero spin-spin cor- 
relations that decay exponentially with distance. These 
properties are shared by other resources previously ap- 
pearing in the literature [H [10j [20l [20] , notably those 
based on the so-called AKLT model [2S [22 ■ We also 
how ID N — U — N chains can be coupled together to 
produce a quasi-deterministic 2D resource for universal 
MBQC. 

Next, we show that when B-type operators act on 
alternating qubits in a ID cluster state, the result is 
in general a probabilistic resource for single-qubit rota- 
tions. We call these states B — U — B chains. We find 
that under a restricted subset of B-type operators, the 
three-dimensional analogs of B — U — B chains constitute 
quasi-deterministic resources for MBQC under strictly 
projective measurements. A similar result was exhibited 
in Ref. [TS], in which a 2D cluster state deformed by 
B-type operators was shown to be reducible to a perco- 
lated 2D cluster state [23] by the action of three-element 
POVMs. The B — U — B states have the interesting prop- 
erty that each qubit can be arbitrarily locally pure, or 
alternatively that an individual qubit can be arbitrarily 
weakly entangled with the rest of the state, as measured 
by the von Neumann entropy of entanglement. Like the 
cluster states, they also exhibit vanishing long-range cor- 
relations, with no spin-spin correlation functions beyond 
second-nearest-neighbor surviving. 

The structure of the paper is as follows: in Section [TTJ 
we briefly review the theory of measurement-based quan- 
tum computation using cluster states, and introduce the 
various definitions and notation used in the technical part 
of this paper. In Section |III[ we describe the effects of 
invertible local operators acting on a cluster state on the 
class of linear transformations that can be logically imple- 
mented via adaptive single-qubit measurements on this 
new state, and outline some strategies for dealing with 
these effects. In Section |TV[ we provide explicit examples 
of some structures of SLOCC-transformed cluster states 
that are universal for either probabilistic or determinis- 
tic single-qubit rotations or full MBQC. Finally, in Sec- 
tion [Vj we discuss the relationship of our resource states 



with previously known quasi-deterministic resources, and 
outline the prospects of identifying hitherto unknown re- 
source states by this method. 

II. BACKGROUND 

An n-qubit cluster state can be defined in terms of the 
stabilizer formalism |24j as the unique n-qubit pure state 
|C1„) satisfying the n conditions 

X, (g) Z,-|CU = |Cl n ), (2) 

where i £ {1, ...,n} labels a qubit, j\f(i) denotes the 
spatial neighbourhood of qubit i, and Xj and denote 
the standard single-qubit Pauli operators, given in the 
computational basis by 



acting on qubit i. Alternatively, the cluster state can be 
identified as the result of a dynamical process in which 

1. n qubits are initialized in the state |+)® n , where 
|+) = (|0> + |1)) is the +l-eigenstate of X; 

2. CZ entangling gates, whose action on the compu- 
tational basis states is given by 

GZ^jlx, y)ij = (-l) x ' v \x,y)ij, (3) 

are applied between each pair of neighbouring 
qubits Thus, 

|C1„) = JJ CZ w |+) l8B , (4) 

where indicates that i and j label neighbour- 
ing qubits. 

For notational convenience, define a global entangling 
operation on a lattice, 

z-i 

®k,f. II (5) 

J= k 

to be the tensor product of CZ gates acting between 
all nearest-neighbour pairs of vertices on a line with la- 
bels between k and I. Now consider a modified one- 
dimensional n-qubit cluster state, 

\c^ D y = (s ljn \m+)t~n- ( 6 ) 

where the first qubit was encoded in some general pure 
state \ip) before the global entangling operation ©i in was 
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performed. The effect of projectively measuring the first 
qubit, the one on which was initially encoded, is to 
teleport the quantum information corresponding to the 
state \tp) to the next qubit, subject to some linear trans- 
formation depending upon the basis and outcome of the 
measurement. To see this, assume that \ip) — a\0) +b\l). 
We then find that 

|CO' = ® 2 „ l CZ 1 , 2 ( a |0+) 1 , 2 +6|l+) 1 , 2 )|+)^ : „ 2 
= «2,«(a|0+) lj3 + 6|l-)i, 2 )|+>^ 



= e 2) n(a|0)iI 2 + 6|l)iZ 2 )|+)S 



(7) 



Projecting the first qubit via an arbitrary rank-1 pro- 
jector |m)(ra|, the state of the system (neglecting the 
projected qubit and overall normalization) becomes 



m = © 2 ,„ (a(m\0)h + b(m\l)Z 2 ) |+)f«^ 

= C5 2 ,„ (a(m\0)\+) 2 + &<m|l>|->a) \+)tJ 

= 2 ,„H 2 (a(m|0)|0) 2 + 6<m|l>|l> 2 ) |+)|^ n 2 

oc 2 ,„H 2 «m|+)I 2 + (m|-)Z 2 ) |^) 2 |+}f;; n 2 , 

where the Hadamard operator Hj = (Xj + Zj)/v2- In 
other words, the quantum information has been tele- 
ported to the second qubit through the linear transfor- 
mation 



M = H((m|+)J- 



\-)Z) 



(8) 



Without loss of generality, the single-qubit state acting 
as the projector can be written as 



/Mi."); =cos||+)+e^sini|-) 



(9) 



for some < £ < 2ir, —it < <fi < %. Thus, the linear 
transformation through which the quantum information 
is teleported can be written as 



\\ II I cos-I + e i0 sin-Z 
2 2 



(10) 



In the special case that 4> — ±7r, corresponding to |m) 
lying on the x — y plane of the Bloch sphere, this transfor- 
mation becomes the (familiar from cluster state MBQC) 
unitary transformation HR 2 [±£]. Thus, there is an entire 
single-parameter family of unitary gates through which 
the initial state \ip) can be teleported, each corresponding 
to a projection of the first qubit on to some state lying 
in the x — y plane. This family is universal for single- 
qubit rotations: via four projections, corresponding to 
£ = 0,£ 2 ,£3,£4 respectively, one teleports the transfor- 
mation 

U(£ 2 ,£ 3 ,£ 4 ) = HR Z [&] HR Z [&] HR Z [6] HR 2 [0] 

= R* [&] [6] B« [6] , (11) 

which is an arbitrary single-qubit unitary decomposed in 
terms of Euler angles. 



To this point, we have not discussed how to com- 
pensate for the randomness associated with the mea- 
surements. If one were to drive the gate teleportation 
described above via projective measurements, measure- 
ments must be made in an orthonormal basis containing 
|m(£, 4> = 7r)). For a single qubit, this basis would be 
B (£, 4>) = {\m (£, 0)) , \m x (£, 0)) } , where 

|m(C,0)> = co S |l+>+e < *sin||->; (12) 

|m x (£,0)) = 8 in||+)-e^cos||-). (13) 

From Eq. (JsJ) , the teleported gate associated with the 
application of the projector |m. )(m | on the first qubit 
would be 



= H sin fl - cos fz 



: XTI cos -I - e-^sin-Z 
2 2 



where in the last step we have made use of the identity 
XH = HZ, and have dropped an unimportant overall 
phase. Once again considering the special case that (f> = 
±7r, this reduces to 

M 1 - = XHR Z [±£] = XM. 

The teleported gate can be summarized succinctly as 
follows: measuring in the basis {Im),!?^)} defined by 
Eqs. ( 12|13 l with <f) — 7T , an d denoting the measurement 
outcome by m — for state |m) and m = 1 for |m ), 
then the teleported gate is X m HR z [£]. The operator X 
can be thought of as a byproduct operator that occurs 
as a result of obtaining measurement outcome 1. 

If the aim is to teleport the operator U (£ 2 ,£3,£/t) de- 
fined in Eq. (Ill, then apparently one runs into a problem 



should a measurement outcome of 1 be obtained for any 
of the four measurements needed to teleport this gate. In 
fact the X byproduct operators can be pushed through 
the rotations because R z [£] X = XR 2 [— £]. Suppose then 
that one performs four projective measurements on a one- 
dimensional cluster state, with the ith measurement be- 
ing a projective measurement of qubit i in the orthonor- 
mal basis B (0j, 7r), 6\ — and the measurement outcome 
denoted rrij € {0, 1}. The quantum information orig- 
inally situated on qubit 1 before the global entangling 
operation is then teleported through the gate 

M = X m4 HR z [0 4 ] X" l3 HR z [0 3 ] X™ 2 HR Z [0 2 ] X mi H 



= x m4 Z m3 X m2 Z mi R T 



(-1) 



1113+mi 



x R 2 [(-l) m2 3 ]R,[(-l) mi 2 ] 



(14) 



Comparing Eq. (11) and Eq. (14 1, the choices 2 = 

£4 make 



(-1)" 11 6, 03 - (-I)™ 2 6, and 4 = (-lH+»n 



the implemented teleported gate 

M = X m4 Z" l3 X m2 Z mi U (£ 2 , £3, £4) 



(15) 
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Thus, any gate can be implemented by conditioning each 
of the last three measurement bases on the results of pre- 
vious measurements, up to an overall Pauli byproduct 
operation. The byproduct is of no concern, as Z has no 
effect on computational basis states while X merely swaps 
them; this means that the effect of the byproduct can be 
taken into account simply by appropriate reinterpreta- 
tion of the final measurement outcomes of the circuit, 
contingent on the intermediate measurement outcomes. 

An equivalent description of this universal gate telepor- 
tation can be obtained within the Matrix-Product State 
(MPS) representation [H [5] of the one-dimensional clus- 
ter state: 

|C1„) = ]T A™ [iJAM fe] . . . AM [i n ] \i t i 2 . . . i n ), (16) 

i 

where i is an n-bit string and the site matrices {A^ [i,]} 
are all two-by-two, except for the boundaries; the 
are row vectors and the {AM[i n ]j are col- 
umn vectors. The site matrices are not unique, but it 
is particularly convenient if they are chosen to satisfy 
the relation E t ^ [iI fe]^ ljlt fe] = 1 for cach h 

corre- 
sponding to the canonical form of the MPS [35] . For the 
left and right boundaries one obtains A^[0] = 4s (+|, 

AW[1] = ^(-\, AW[0] = |0), and A^[l] = |1); for 
the bulk sites 1 < j < n they are A^[0] = ^75 H and 
AW[1] = ^HZ = ^XH. The 'always-on' operator H 
is teleported on each measurement of a qubit, and the 
X gate serves as the byproduct operator. In general, an 
MPS state is a universal resource for measurement-based 
single-qubit gate teleportation (a 'computational wire') 
if the bulk site matrices can be chosen to be propor- 
tional to unitaries [TS]. In this case they can be writ- 
ten as AW[0] = -^W and AW[1] = -^WR z (cj)) with 
W <= SU(2) and 4> <= E. 

In the context of the calculations presented in the 
next two sections, it is worth pointing out that there are 
two special features associated with the projective single- 
qubit measurements on one-dimensional cluster states 
presented above. The first is that the linear transfor- 
mation M on the quantum state \tp) is guaranteed to be 
unitary; in practice, this means that the effect of such 
a measurement is not dependent on the input state \if)). 
If the linear transformation is not unitary (for example 
projections of the local system outside the x — y plane, 
as discussed below), then there is an equivalent unitary 
transformation resulting in the same final state vector; 
however, the equivalent unitary will depend on \ip). The 
second feature is that the byproduct operator resulting 
from measurement outcome 1 is always X. This is bene- 
ficial as X operators can be pushed through R z operators 
with an easily characterized effect, as discussed above. 
These properties can be summarized as follows: 

• Property IA: the teleported gate is in general of 
the form HR Z [£] where i.e. a unitary gate 



corresponding to a z-axis rotation by a real angle, 
followed by a Hadamard gate. 

• Property IIA: the byproduct operator is always 

x = r, [4 

The above two features do not hold for single-qubit 
projective measurements outside the x — y plane. In gen- 
eral, 

• Property IB: the teleported gate is in general of 
the form HR Z [£] where ( € C, i.e. a non-unitary 
gate corresponding to a z-axis rotation by a com- 
plex angle, followed by a Hadamard gate. 

• Property IIB: the byproduct operator is in gen- 
eral K x [rj] where 77 7^ tt in general, i.e. an x-axis 
rotation not corresponding simply to Pauli X. 

Another way to view this is that in x — y plane, one 
always teleports HR Z [£ + S mt ±%] with £ G K, whereas 
in any other plane, one teleports HR Z [£ + <5 m ,ie] where 
£ € R corresponds to the angle of the z-rotation in the 
desired gate, and e € C is a complex error that occurs on 
measurement outcome 1. 

Single-qubit gates alone are not sufficient for universal 
computation; at least one multiqubit entangling gate is 
required as well. In the cluster state model, multiqubit 
gates are accomplished via measurement patterns on 2D 
structures. Logical qubits are processed by horizontal 
ID wires, while entangling operations are mediated by 
vertical links between them. An entangling gate that is 
locally equivalent to controlled-NOT can be achieved by 
measuring a link qubit in the Y basis [IJ [5] . 

As with the case of single-qubit rotations, local Pauli 
byproduct operators may exist as well, depending on the 
measurement outcomes. As before, these byproduct op- 
erators are of no concern computationally. Thus there ex- 
ists for cluster states a measurement pattern that deter- 
ministically implements a two-qubit unitary entangling 
gate, with any byproducts that occur being of the local 
Pauli type. It is not immediately clear that this will be 
the case for states other than the cluster state; in gen- 
eral, the teleported two-qubit gate may be non-unitary 
and the byproduct may be non-local. Any candidate re- 
source state for MBQC must be shown to be amenable to 
a measurement pattern implementing some suitable two- 
qubit entangling gate. In Section |III[ we describe some 
ID structures that are resources for single-qubit rotations 
and then in the examples of Section |Tv] we demonstrate 
how to perform entangling gates with natural 2D or 3D 
extensions of the ID structures, and how to compensate 
the randomness associated the measurements. 



III. PROJECTIVE MEASUREMENTS ON 
SLOCC-TRANSFORMED CLUSTER STATES 

As discussed in the previous section, the distinguishing 
feature of MBQC with regular cluster states is that there 
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exists a plane of the Bloch sphere onto which successive, 
adaptive, single-qubit projective measurements drive an 
arbitrary computation that is deterministic and of fixed 
length. No matter which single-qubit gate is desired, it 
will be implemented with certainty up to an unimpor- 
tant Pauli byproduct with four measurements. For an 
SLOCC-transformed cluster state, it is not obvious that 
there exists any such plane: in general it is not possible to 
simultaneously satisfy both Properties IA and IIA (or for 
the latter any another convenient Clifford gate). A nat- 
ural question to ask is then: under what circumstances 
can either property IA or IIA be satisfied by itself? And 
if only one property is satisfied, does there remain a de- 
terministic protocol for universal quantum computation? 
Sec. MI Al and Sec. IIII 51 discuss the circumstances under 
which it is possible to independently satisfy Property IA 
and IIA, respectively. 



or equivalently 

|m) = v^(5 (1)t ) _1 e~ ia 

= e- ia (sWy 1 R z [$\+) 



cos 1 1 



i sin - 
2 1 



(21) 



Eq. (21 1 is the condition on the state \m) such that 
measurement outcome yields a unitary teleported gate. 
Note that there is a family of states characterized by a 
single parameter £ fulfilling this condition, not including 
the unimportant overall phase a. 

To ensure that the measurement yields a unitary tele- 
ported gate independent of the measurement outcome, a 
similar condition must follow for the orthogonal comple- 
ment | m ). Orthogonality requires 



(m x \ot(-\R z [-£] fiWt, 



(22) 



A. Strategy I: Guaranteed Unitary Evolution 

1. Derivation of N -type Operators 
For convenience, define 



Sfc.z 



(i) 



(17) 



j=k 



where G GL(2,C). From Eq. Q, it is clear that 
the SLOCC-transformed cluster state encoding quantum 
information can be written in the form 



i I | \<g>n— 1 
2 ) l+/2,...,n- 



(18) 

Following the procedure discussed in Sec. [TT} applying the 
projector |to)(to| to the first qubit yields the resulting 
state on the remaining qubits: 



S 2 ,„© 2 ,„ (a{m\s[ l) \0)l 2 + &<m|£f 
S 2 ,„© 2 ,„H 2 M 2 |V>) 2 |+)!rn> 



\ (girt— 1 
12... .Jn 



where 

M 2 =H 2 



v/2 



-I 2 + 



i\SW\ 



V2 



(19) 



The only way for this to correspond to a unitary gate is 
if 



^Hs (1) h> = - 

where < £ < 2tt, and therefore 
S (1)t |m) = V2e- ia 



— e cos 



-ie la sin|, 



£ £ 
cos-|+) +isin-|-) 



(20) 



and therefore 



cS {1) R z [£] 



(23) 



for some c ons tant c G C. Repeating the procedure that 
led to Eq. (20), but with \m^-) instead of |to), one obtains 



£ 



ism -| 



(24) 



Substituting Eq. (|23|) into Eq. (|24|) yields 

= cS^sWR, [£] 



cos ■ 



£ 



■ism-\-} 



(25) 



Rewriting Eq. ( 25 ) in the computational basis results in 
the expression 

|0) + e-*|l) = c'S (1)t S (1) (|0) - e 4 «|l)) (26) 

for a suitably defined constant c'. Defining := 
g(W S (k) and T i,j ._ (i\T k \j), we can see from Eq. Q 
that 



1; 



From here, it is easily deduced that 



T o,o _ gl 



0.1 



1.0 



(27) 
(28) 

(29) 



Making use of the Hermiticity of simple algebra yields 
T o,o =Tj u (go) 



Eq. ( 30 1 above has a very simple geometric interpreta- 
tion: it means that must preserve the relative norm 
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of the computational basis states. There is no require- 
ment for S^ 1 ' to preserve their orthogonality, however, 
which means that is allowed to differ quite drasti- 
cally from a unitary transformation; in fact, it can be 
made arbitrarily close to singular, as the relative angle 
between the computational basis states under the trans- 
formation by S^ 1 ' can be vanishingly small. We will refer 
to operators obeying this norm-preservation restriction 
as N-type operators. 

Definition III.l. A GL(2,C) operator S satisfying 
(0\S'S\0) = (1\S* S\l) is called an N-type operator. 

The singular value decomposition is helpful for charac- 
terizing N-type operators. An arbitrary SLOCC operator 
S can be written in terms of its singular value decompo- 
sition as S = U-DV, where U is an arbitrary two-qubit 
unitary, D is a positive-definite diagonal matrix 



and are found to be 



D = K 



cos 9 
sin( 



(31) 



where < 9 < ~, and V is an arbitrary unitary ma- 
trix parametrized via the Euler decomposition as V = 
HR 2 [a] R x [P] R z [7], < a, /3,7 < 2n, and any global 
phase has been absorbed into U. It is then straightfor- 
ward to determine that 

(0|S' t S'|0) = ^k 2 (1 + cos 26> sin a sin/?); (32) 



(l|StS|l) 



- k 2 (1 — cos 26* sin ot sin j3 ) . (33) 



So, if S is an N-type operator, we must have 9 — ? (in 
which case S is proportional to a unitary), a = or 
f3 = 0. The case where f3 = still allows us to assume 
a = without loss of generality. Doing so, the [/?] 
operator can be commuted past the H to turn into a z- 
rotation, and then absorbed into U. Thus we have the 
following characterization of N-type operators. 

Lemma III. 2. Every N-type operator S must either 
be proportional to a unitary operator, or of the form 
S = UDV , where U is an arbitrary two-by-two unitary 
operator, D is defined as in Eq. (31) and V = HR 2 [7] 
with < 7 < 2tt. 

It is straightforward to obtain an expression for the 
byproduct angle // in the case of measurement outcome 
1 in terms of the parameters 9, 7 and £ (recall that this is 
the degree of freedom in the measurement basis). Using 
Eq. (19), it can easily be checked that the teleported 
gate when \m) oc (S^) -1 R z [—£] H|0) is M = HR Z [£], 
and when \m) oc SR Z [-£] Hjl) is M 1 - = R x [fjf] HR 2 [£], 
where the byproduct angle \j! obeys 



u! 1 — cos 26* cos (7 — E) 
tan ^- = 



cos 29 sin (7 — £) 



(34) 



The probabilities of the two measurement outcomes can 
also be easily calculated in terms of the same parameters, 



P(0) = 2 (1 + cos 26» cos 20 



P(l) 



(I - cos 26 cos 20 



(35) 
(36) 



This differs from the case of gate teleportation with a 
perfect cluster state, where the two measurement prob- 
abilities are always exactly |. That said, the expected 
probability of obtaining a byproduct here, averaged over 
all £, is 



(p(i))e = \, 



(37) 



irrespective of 9. We can thus generically expect to ob- 
tain an unwanted byproduct operator that we must com- 
pensate on half of our single-qubit measurements. This 
point will be discussed in Sec. |IV| 



2. Properties of N- Transformed Cluster States 

Cluster states locally transformed by N-type operators 
can exhibit remarkably different properties from perfect 
cluster states. Nevertheless, as will be shown later in 
Sec. IV B[ they can under some circumstances serve as 
universal resources for MBQC of random length. 

Consider for the moment the Schmidt decomposition 
of an n-qubit cluster state on some set of qubits V with 
respect to a bipartition separating qubit k from the rest: 



ICln), 



1 

75 



(|0)fc|Cln_i) v \fc + |l)fcZ iA /-( fc )|Cl n _i)v\fc) 



(38) 

where |Cl„_i)y\fc refers to the cluster state resulting from 
deleting qubit k and Z^( fc ) is the tensor product of Z op- 
erators acting on all the neighbors of k. This can be 
checked by verifying that this state satisfies the stabi- 
lizer conditions ([2]) . The Schmidt basis for the multiqubit 
component can be further decomposed if desired by the 
same technique. The equality of the Schmidt coefficients 
in Eq. (38 1 demonstrates that any individual qubit in a 



cluster state has a maximally mixed local reduced density 
matrix, or in other words that it is maximally entangled 
with the rest of the cluster with respect to the von Neu- 
mann entanglement entropy. Likewise, exactly one ebit 
of entanglement is shared across any bipartition of the 
cluster state. 

One effect of N-type operators is to change the lo- 
cal reduced density matrices of individual qubits within 
the state. In the canonical representation, the site 



matrices of the MPS representation [cf. Eq. ( 16 )] are 
A0[O] = -^W = HR Z [— (7O) +20W)] and I%] = 

^7jIUR z [40")], which are both unitary. This immedi- 
ately implies that the channel having these matrices as 
Kraus operators is unital, and like the ordinary cluster 
state one ebit of entanglement is shared across any bipar- 
tition. That said, the entanglement between any given 
qubit and the rest of the system need not be unity. 
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Consider for example the local reduced density matrix 
of a qubit adjacent to an endpoint of a ID cluster state 
with n qubits, numbered 1 to n from left to right. Sin- 
gling out first qubit 2 and then qubit 3, this state can be 
written as 



I CI 



1D\ 



1 



v/2 



(|0)s 



ID 



>iZ 3 



|C&3. 



(39) 



Now consider the action of an N-type operator, — 
DHR Z [7] on qubit 2, where the leading U operator is 
dropped because it can be absorbed into the measure- 
ment basis. It is easy to check that 



V2 



(cos0|O) 2 |$)-isin0|l) 2 |$ x )) 



where |<&) and j^- 1 ) 
1,3, ... ,n such that 
a Schmidt decomposition. The local reduced density ma- 
trix of qubit 2 is 



(40) 

are 7-dependent state s fo r qubits 
= 0. Thus, Eq. po} remains 



(>2 



cos 2 9 








revealing that qubit 2 is no longer maximally entangled 
with the rest of the state. A similar calculation can be 
performed for qubits further from the boundary, with 
qualitatively similar results. 

Another property of these states is the long-range be- 
havior of two-point correlation functions, those of the 
form Cij(A, B) := (AiBj) — (Ai)(Bj) for some oper- 
ators A and B. Two-point correlation functions of 
large ID cluster states with periodic boundary condi- 
tions can be efficiently calculated using the Matrix Prod- 
uct State (MPS) representation. For an n-qubit ring, 
calculation of the correlation functions amounts to tak- 
ing traces of products of n 4 x 4-dimensional matri- 
ces. Consider therefore a ID cluster state with peri- 
odic boundary conditions (i.e. a ring), with the opera- 
tion 

N (i) 

acting on qubit i. For this state, calculations 
show that all two-point Pauli correlation functions van- 
ish except for the second-nearest-neighbour correlation 
function Ci_i,i + i(Z, Z) = cos2#sin7. This is in con- 
trast to the perfect cluster state with periodic boundary 
conditions, for which all two-point correlation functions 
identically vanish. 

As another example, the relevance of which will be- 
come clear in Sec. IV A| consider a ring with an even 
number of qubits, with N acting on every alternate qubit; 
say, the ones with even labels. In this case, the magni- 
tude of the same two-point correlation function between 
odd-numbered qubits decays exponentially: 



\C 1>2j+1 [Z, Z)|^ exp (-^ 



(41) 



The length scale L depends on 9 and 7. The same corre- 
lation function between pairs of qubits with at least one 
even label is zero. The numerically obtained behavior of 







71 












4n 


L 


3 n 


2.5 




^ 40 












K 


2 




8 






In 












40 


1.5 




9ir 








1 




40 


0.5 








FIG. 1: (Color online) Correlation length scale L associated 
with the correlation function \C\,2j+i{Z, Z)\ ~ exp (— for 
a N — U — N ring with all N-type operators identical, on a ring 
of 1000 qubits, as a function of the parameters of 7 and 8, the 
parameters of N. The length scale increases as N approaches 
the singular limit, i.e. as 6 gets close to or 7r/2. 



L for a ring of 1000 qubits is shown in Fig. [T] as a func- 
tion of 7 for several values of 9 between and 7r/4. As 
can be seen from the figure, the length scale increases 
with decreasing 9 over this range, i.e. as the N-type op- 
erators approach the singular limit 9 = 0. The length 
scale is symmetric about 9 = 7r/4 between and ir/2. 
Viewed as a function of 7 with 9 held constant, the cor- 
relation function is convex and non-negative in 7 over 
the interval from to tt/2 and is symmetric about tt/4. 
For 7r/2 < 7 < 7r, the magnitudes behave the same way 
as in the previous interval, but the signs alternate. The 
7-behavior is periodic with period it. We note in passing 
that these non-zero correlation functions provides a lower 
bound for the localizable entanglement [26] between that 
pair of qubits in the state via projective measurements, 
with respect to the concurrence [27j . 

Note that a number of resources for MBQC with 
non- vanishing long-range correlation functions have been 
pointed out in the literature [TUl HI HI 120], based 
on the so-called spin-1 AKLT model [22 E2J- These 
states are quasi-deterministic resources, in the sense that 
measurement-based computations using these states can 
be made arbitrarily likely to succeed, either by reduc- 
tion of the resource state to a deterministic resource or 
by a repeat-until-success strategy with each elementary 
gate requiring a random number of measurements. In 
Sees. |IV A| and |IVB[ we describe resource states called 
N — U — N states that are based on cluster states trans- 
formed by N-type operators; these states share the prop- 
erties of quasi-determinism and non- vanishing long-range 
correlations. 



B. Strategy II: Guaranteed Pauli Byproduct 

1. Derivation of B- Type Operators 

Another possible strategy is to attempt to ensure that 
the byproduct operator is guaranteed to be Pauli-X, 
whether or not the tclcportcd linear transformation is 
unitary. The advantage of this approach is that X has 
nice commutation properties through rotation operators 
about the z-axis, whether they be by real or complex an- 
gles, leading to the hope that the randomness inherent 
in the measurement process can be easily compensated. 

When projectively measuring in the orthonormal basis 
{|m), jm" 1 "}}, the two possible operations that can occur 
are 



M = H 
M x = H 
= XH 



i\SW\ 



V2 
i x \S^ 



-I 



m x \S^\ 



V2 V2 

V2 V% 



(42) 
(43) 
(44) 



Since I and Z are linearly independent, it follows that 
for the byproduct to be guaranteed to be proportional to 
Pauli-X, one must have 



or equivalently, 



(m\S\+) 
(m\S\-) 



(m\S\0) 
(m\S\l) 



c(m x \s\-y, 

c(m x \S\+), 

c(m- L |5|0); 
-c{m x \S\l), 



(45) 
(46) 



(47) 
(48) 



for some non-zero constant c G C. Suppose S = U-DV 
where U is an arbitrary single-qubit unitary, D is de- 
fined in Eq. Q, and V = R z [(3} R x [7] R z [S\. Further 
UU'\0) and |m x ) = UU'\1), with 
The reason for the appear- 



suppose that \m) 
U> = R z [13'} R x [7'] R z [S'} 
ance of U in the definitions of |m) and \m x ) is to com- 
pensate for the appearance of U in the singular value 
decomposition of S. The only effect of the R z [8'} opera- 
tion is to multiply the teleported gate by a global phase, 
so we can choose 8' = in what follows without loss of 
generality (it remains a free parameter for the applied 
unitary U'). Having done so, Eqs. ( 47p8 | can be rewrit- 
ten as 



(0|Q|0) = c(l|Q|0); 



(0|Q|1) 
where we have defined 
Q := (U'^DV 



-c(l|Q|l) 



Rx [-7'] Rz [-13'} DR Z [(3} R x [7] R z [8] 
R x [-7'] Rz [f3 - 13'} DR X [7] R z [8} 
Rx [-7'] Rz [b] DR X [7] R z [8} ■ 



(49) 
(50) 



(51) 



(52) 



in the last line above we have defined b := /3 — (3'. In the 
expression above, 7' and b are free parameters, while D, 
7 and 8 are determined by the SLOCC operator S. 

Return now to the constraints, Eqs, ( 49|50 ). Denoting 
Qij ■= (i\Q\j), one finds that Q00/Q10 = -Qoi/Qn- 
Note that neither Q10 nor Qn can be zero; if either were 
zero, then the constraints would force Q and therefore S 
to be singular, which by assumption is not the case. This 
in turn means that 



Det(Q) = 2QooQii. 



(53) 



From the definition of Q, Eq. (52 I, 



Det(Q) = sin 2(9; 
2QooQn = sin 7' sin 7 (cos 7 — i cos 26 sin b) 

+ sin20 (1 + cos 7 cos 7') . (54) 



Substituting the above expressions into Eq. ( 53 ) and 



equating real and imaginary parts gives us the two con- 
ditions 

sin 7' sin 7 sin b cos 29 — 0; (55) 
cos 7 (sin 7' sin 7 + cos 7' sin 26) = 0. (56) 

Recall that since S is invertible, we cannot have sin 26 — 
and since S is non-unitary, we cannot have cos 26 = 0. 
The only ways to satisfy Eq. ( |55"|) ar e if sin 7' sin 7 = or 
sin 6 = 0. In the first case, Eq. ( |56[ ) immediately implies 
that cos 7' cos 7 = 0, leaving b as a free parameter for 
our measurement basis. In the second case, 7' is fixed in 
terms of 6 and 7, leaving no freedom in the measurement 
basis we are using. Furthermore, if we choose sin 7' = 0, 
i.e. 7' € {0, 7r}, then the measurement basis we are using 
is restricted to being the computational basis acted on by 
U (completely specified by S) ; again, no freedom. There- 
fore, the only solutions available to us that leave freedom 
in the measurement basis, and thus the teleported gate, 
are 7 € {0,tt} and 7' G {§ , ^}. Note that 



5 f 5 = 



1 + cos 7 cos 26 



sin 7 cos 26 



' sin 7 cos 26 1 — cos 7 cos 26 



(57) 



Thus, demanding that the SLOCC operators allow for a 
guaranteed Pauli by-product, assuming the SLOCC op- 
erator is not unitary and thus cos 26 ^ 0, is equivalent 
to demanding that S'S be diagonal. Geometrically, this 
means that S must preserve the overlap of the computa- 
tional basis states (the transformed computational basis 
is still orthogonal). We will refer to this kind of basis- 
preserving operator as a B-type operator. 

Definition III. 3. A GL (2, C) operator S satisfying 
(0|S' t S'|l) = (llS^SlO) = is called a B-type operator. 

A B-type operator can therefore be written 



B = 



XJDR z [p}R z [8}, 7 
UDR z [j3}XR z [8], 7 
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ignoring overall phases. The two possibilities above can 
be simplified and collapsed into one. First, note that 
R z 1(3] can be commuted past D and absorbed into U. 
Next, note that XDX is itself a diagonal matrix that 
results from swapping the diagonal entries of D. This 
means that the case where 7 = ir can be written instead 
as U'D'R 3 [5], where U' = UR Z [(3] X and D' = XDX. Of 
course, R 2 [S] can also be absorbed into U'; thus, a simple 
and completely general expression for a B-type operator 



B = UL>. 



(58) 



The diagonal matrix D in the singular value decomposi- 
tion can be expressed as 

D cx diag(cos(0),sin(0)) = y/sm{0) cos(0)R 2 [i In cot(0)] , 
so that the B-type operator becomes 



B cx UR z [ilncot(0)]. 



(59) 



Because the unitary U can be absorbed directly into the 
measurement basis, one can interpret B-type operators 
as z-rotations by an imaginary angle, the value of which 
is related to the ratio of the singular values. 

When the local operator is B-type, the single- 
parameter family of measurement bases satisfies 7' € 
{§>!?}> an< ^ P' e [0, 27r) is a free parameter. When 
this family of bases is used, the byproduct operator as- 
sociated with measurement outcome 1 is always Z (up to 
a global phase). The teleported linear transformation is 
no longer unitary, however; it takes the form of a rota- 
tion about the z-axis of the Bloch sphere by a complex 
angle, followed by a Hadamard operation. The real part 
of the angle is completely specified by the choice of mea- 
surement basis, via the free parameter (3 1 . The imaginary 
part is purely a function of the ratio of the singular values 
of the local GL(2,C) operator. Denoting the measure- 
ment outcome corresponding to 7' = ? by m = and 
that for 7' = y by m = 1, the teleported gate is given 
(up to a global phase) by 



M = X m HR 2 [/?' + iln cot 0] 



2. Properties of B-Transformed Cluster States 



(60) 



Interpreting the B-type operators as z-rotations by 
imaginary angles provides a simple insight into the nature 
of B-transformed cluster states. The R 2 operator com- 
mutes with all CZ gates, so one can push it all the way 
through to the |+) states in the definition of the cluster 
state, Eq. Q. Because R 2 (£)l+) is an arbitrary single- 
qubit state, B-transformed cluster states are equivalent 
to applying CZ gates between qubits in arbitrary states 
(not including computational basis states, which would 
require singular B operators). 



One might assume that B-transformed cluster states 
are equivalent to weighted cluster states [28H32] , but this 
is not in fact the case. Weighted graph states are defined 
as ri(i,j) CP(</?)ij 1+)®", where the controlled-phase en- 
tangling gate is CP (ip) — diag (l, 1, l,e lv ); the cluster- 
state edge weights are then given by aiy = ipij. Con- 
sider the simplest counter-example of a three-qubit lin- 
ear cluster state with the central qubit transformed by a 
B-type operator B = DR Z [7] with D — diag(cos0,sin0). 
The eigenvalues of the local reduced density matrices are 
all {h (1 ± cos 20)}. On the other hand, for a three- 
qubit ID weighted graph state with edge weights (fi2 
and y>23, the eigenvalues of the reduced density ma- 
trix are ^(licos^ 1 ), | (l ± | cos <fx2 cos ^23) , and 
\ (licos^ 2 ) for qubits 1 through 3, respectively. If 
the weighted graph and the B-transformed cluster are 
LU-equivalent, there must be some choice of ifi2 and 
if '23 such that the spectra of the reduced density ma- 
trices are the same in both cases. For qubits 1 and 
3 this implies <j>i2 = <p23 = 40. For qubit 2 one ob- 
tains I (l ± I cos40 2 ). This matches the correspond- 
ing spectrum for the B-transformed cluster only when 
= ±j,tp = ±7r, in which case both states are LU- 
equivalent to a perfect cluster. 

Cluster states locally transformed by B-type opera- 
tors also exhibit different properties from perfect clus- 
ter states. As with N-type operators, B-type operators 
change the local reduced density matrices of individual 
qubits within the state, as described in the following 
lemma. 

Lemma III. 4. Let |Cl n )v be an n-qubit cluster state 
on the set of qubits V, with some subset Q C V acted 
upon by P-type operators. In particular, suppose that for 
each qubit i G Q, the P-type operator acting is given by 
Ed) = £>(<) with £>« = \/2diag (cos 0^, sin 0^). Then, 
the local reduced density matrix for any qubit k G V is 
given by 

Pk = cos 2 (fc) |O)(O|+sin 2 0( fc )|l)(l| 
sin 

20{k) 11 cos20W|O}<l|+h.. 

where 0^ := f if k £ Q. 

The lemma is easily proved by taking advantage of the 
expression ( 38 ) for the Schmidt decomposition of a clus- 
ter state with one subsystem being qubit k alone, and 
then calculating p k directly. The calculation is done by 
expressing the cluster state as the action of controlled-Z 
gates acting on the product state |+)® n , and then using 
the fact that the £)W and controlled-Z gates are mutu- 
ally commuting. A consequence of this lemma is that 
the reduced density matrix of a given qubit is maximally 
mixed if and only if the qubit itself and at least one of 
its neighbors are untouched by B-type operators. In gen- 
eral, qubits within B-transformed cluster states are not 
maximally entangled with the rest of the state; in fact, 
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they can be arbitrarily weakly entangled (with respect to 
the von Neumann entanglement entropy). 

Consider now a ring of an even number of qubits, with 
identical B-type operators specified by #( 2fe ) — act- 
ing on the qubits with even labels. The significance of 
such a state will become clear in Example |IVD| Two- 
point correlation functions can be calculated exactly for 
this state. The result is that the nearest-neighbor cor- 
relation functions C2fe,2fc±i(Z, X) = cos 29 and the next- 
nearest-neighbor correlation functions C2fc,2fe±2(Z, Z) = 
C*2A;-i,2fc+i(X, X) = cos 2 (2#) are the only ones that are 
non-zero, while all the other two-point Pauli correlation 
functions are identically zero. Again, this differs from 
the perfect cluster state, where all correlation functions 
are zero. 

Exact calculations on chains of up to 7 qubits 
with identical (but arbitrary) B-type operators B2j = 
DR Z pyW)] acting on even qubits 2j reveal that the 
non-zero Schmidt coefficients corresponding to any bi- 
partition of the chain into two contiguous halves are 
{cos 9, sin 8}. The von Neumann entropy of entangle- 
ment is equal to the Shannon entropy of this list, and 
is generally less than one ebit. The MPS representation 
bears out this observation. In the canonical form the site 
matrices for the boundary qubits are 

-^(cos0( 2 )(O| -t-sin^(l|); (61) 

-^(cos^(0|-si 
AN[0] = |+>; AW[i] =_|_), (63) 



AW[Q] 
AW[1] 



2) (i|); (62) 



those for the bulk even sites are 



A^'][0] = 



e 




%i 





(2j) 







and those for the bulk odd sites are 

cos 6^) sin6»( 2 ^ 



AM-V[0] = 



cos6>( 2j ) sin 9^ 

-cos6»( 2 ^ sin6>( 2 ^ 
cos6K 2 ^ -sin 6^ 



(64) 
(65) 

(66) 
(67) 



It's very easy to verify that the channels induced by the 
matrices on the odd sites are not unital in the sense given 
in Ref. [18] . so a B — U — B chain is not a quantum wire. 

Such a ID state would appear not to be capable of reli- 
ably processing a single qubit. This is true, but a simple 
modification of the geometry from one to two dimensions 
yields a usable resource for random length computation. 



This will be elaborated upon in Example IV D 



IV. RANDOM LENGTH COMPUTATION 

Neither Strategy I nor Strategy II discussed in the pre- 
vious section directly offers a way to perform determin- 



FIG. 2: N — U — N state, a one-dimensional structure that 
can be used for deterministic random-length single-qubit ro- 
tations. 



istic single-qubit rotations. For Strategy I, it is unclear 
how to compensate for a byproduct operator ~R, X [77] where 
rj 7^ 7r, as such a byproduct operator does not possess con- 
venient commutation properties with the H and R z op- 
erations. Similarly, for Strategy II, it is unclear whether 
some number of non-unitary teleported gates can be com- 
bined to form a desired unitary. 

Another perspective on the strategies is that a sin- 
gle measurement with outcome 1 teleports the gate 
HR Z [£ + e] , where £ G C is some angle associated with 
the always-on operation HR 2 [£] (in the terminology of 
Ref. [18] ) and e e C is a possibly complex error associ- 
ated with the byproduct. To correct this error in princi- 
ple requires two additional measurement steps. The first 
measurement step should teleport the gate HR Z [0] = H, 
which would cancel the previously applied Hadamard 
gate; a possible X byproduct operator might result de- 
pending on the measurement outcome. On the sec- 
ond measurement step one would attempt to teleport 
HR Z [— e] or HR 2 [e] depending on the previous measure- 
ment outcome, thus cancelling the original error e. 

This procedure is only possible if the measurement im- 
mediately after first incurring an error cannot itself gen- 
erate any further error e'. One way to guarantee such 
a circumstance is to impose that every alternate Si op- 
erator is in fact unitary. Thus there must exist a class 
of states that are a strict subset of SLOCC-transformed 
cluster states, which constitute resources for random- 
length universal gate teleportation. Likewise, a subset 
of SLOCC-transformed cluster states in two dimensions 
must be universal resources for MBQC. The remainder 
of this section is devoted to various explicit examples. 



A. Deterministic single-qubit rotations: N — U — N 
state 

Consider a one-dimensional state of the form 



\R) = ® U^' 



T (2) ^ N 0) 



uf-.-^N^ICU, (68) 



where the {iV^} are N-type operators, and the -*} 
are local unitaries (c.f. Fig. [2]). The goal is to teleport 
the single-qubit unitary 



u(c,r7,0 = iUC]R* 



(69) 



The first step is to use Strategy I to attempt a telepor- 
tation of HR Z [0], analogously to the scheme with the 
perfect ID cluster state. For the correct choice of basis 
the measurement outcome mj = corresponds to suc- 
cess. One can then immediately measure qubit 2 in a 
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basis that teleports X m2 HR z [£], and then use Strategy 
I to attempt the teleportation of HR 2 [(— l) 1 ™ 2 rj\ starting 
on qubit 3. 

If the first measurement outcome is instead mj = 1 
then one instead teleports HR 2 [e] with e G M. This 
error must be immediately corrected, because the next 
desired rotation is around an orthogonal axis. Happily, 
there is a local unitary U^ 2 ) acting on the next qubit 
in the chain. The Hadamard operator that effects the 
now-undesired transformation of the rotation axes can 
be eliminated by teleporting another one (H 2 = I). This 
is accomplished by measuring the next qubit in the basis 
{lj( 2 )|+), U^ 2 - 1 1— ) } . Labelling the measurement outcome 
777-2 j the teleported gates are 



X™ 2 HHR 2 [e] = X™ 2 R 2 [e] 



(70) 



The measurement basis for qubit 3 is then chosen such 
that measurement outcome 7773 = results in the gate 



HR 2 



(-l) m2+1 e 



being teleported. In this case, the 



overall unitary becomes 



HR Z 



(-l) m2+1 e 



X m *R z [e] = Z m2 HR 2 [-e]R 2 [e] 

= z m2 n. 



At this point one has successfully teleported a Hadamard 
gate and an unimportant Pauli byproduct. The next 
measurement on a qubit with an even label can teleport 
the desired HR 2 [£] gate without error. One then at- 
tempts to teleport HR 2 [77] by measuring qubit 5, using 
Strategy I, etc. 

The procedure corresponds to the following steps: 

1. Measure qubit 1 with outcome mi in the basis 
{(N( 1 )t)- 1 R 2 K 1 ]H|0},N( 1 )R z [-e 1 ]H|l)}; (71) 

2. If mi =0, then success; 

3. If mi = 1 then one has effectively teleported the 
gate R^ [e (1) ] HR 2 [&], where 

m , „ cos 20W C os £1 , m . 

e (1) = ±2arctan s + tt. (72 

l±cos26'( 1 )sina 

Note that eW = when = (6>W = tt/4), 
as expected. Measure qubit 2 with outcome 7772 in 
the basis {U( 2 )X|0),U( 2 'X|1)} ; 

4. Measure qubit 3 with outcome 777.3 i n the basis 
{(N(3)t)-!R 2 [ x ]H|0),N( 3 )R 2 MH|1>}, where X = 

(_ 1 )m 2£ (l). 

5. Repeat steps 3 and 4 on successive qubits 2k and 
2k + 1 until outcome 1 is achieved on an odd qubit, 
using U< 2 ) ->■ \J^ 2k \ N< 3 ) N( 2fe+1 >, 777 2 777 2fc , 

e (l) -> e (2fc-l). 



The key point of this example is that as for any measure- 
ment on an odd-numbered qubit that yields the 'correct' 
outcome m t = 0, one will have succeeded in implement- 
ing part of the desired single-qubit rotation. Further- 
more, any errors resulting from outcomes m, = 1 are 
correctible by making further measurements. This thus 
constitutes a repeat-until-success strategy, and gives rise 
to a quasi-deterministic random-length single-qubit rota- 
tion. The likely reason for this one-dimensional state to 
be capable of processing a logical qubit is that the left 
and right parts of the state share an ebit of entanglement 
with respect to any cut, as mentioned in Sec. |III A 2| 



Deterministic Universal MBQC: 2D N U 
State 



■N 



For universal MBQC, a two-dimensional resource state 
is required. The precise geometry of the two-dimensional 
state on which MBQC occurs is determined by the spe- 
cific circuit to be implemented. Ideally one would start 
with a state defined on a convenient and simple geometry, 
and then 'carve' the desired shape out by deleting certain 
qubits. For cluster-state MBQC, for example, one carves 
the required state out of a rectangular lattice by projec- 
tively measuring the unwanted qubits in the computa- 
tional basis. The goal is to yield isolated one-dimensional 
wires, each of which represents a logical qubit, with links 
only existing between wires in places where an entangling 
gate between logical qubits is needed. 

Consider now a regular two-dimensional lattice com- 
posed of N — U — N states, as depicted in Fig. [3] As in 
the usual cluster state, logical qubits are processed by 
alternating horizontal wires composed of physical qubits, 
and entangling gates by vertical chains connecting them. 
Unlike the cluster case, however, the procedure for im- 
plementing single-qubit rotations with N — U — N states 
is of random length, so it is impossible to decide in ad- 
vance where the desired links between wires will occur. 
The computational cluster state then must be carved 'on 
the fly' 

Suppose that the quantum information encoding two 
logical qubits resides on (yet unmeasured) N-transformed 
physical qubits on two different wires. If an entangling 
gate between logical qubits is not desired at the next step, 
then the link between the wires can be first severed by 
measuring the intervening U-transformed chain qubit in 
the computational basis. An example of the method to 
decouple qubits 1 and 3 is shown in Fig. [3j where qubit 
2 is measured in the computational basis. This has the 
effect of teleporting a Z gate to each of the logical qubits if 
the measurement outcome is m = 1. Other than taking 
into account the possible existence of these byproduct 
operators, the computation subsequently proceeds as in 
the one-dimensional case discussed above. 

The desired entangling gate is implemented as follows. 
At the time that an entangling gate is needed, the local 
part of the resource state looks like two ID N — U — 
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FIG. 3: 2D N — U — N state, universal for quasi-deterministic 
MBQC. The qubits labeled 1, 2 and 3 can be used to imple- 
ment an entangling gate, which involves measuring qubit 2 in 
the Y basis. Alternatively, if an entangling gate is not desired 
here, qubits 1 and 3 can be decoupled by measuring qubit 2 
in the Z basis. 



N states, each coupled via CZ operations to an ancilla 
initially in the state |+) and subsequently acted on by an 
arbitrary U. In Fig. [3j the entangling link is represented 
by the vertical N— U — N chain labeled by qubits 1, 2, and 
3. The local part of the state is mathematically described 
as 



\R) = nPu£ 2) n£ 3) CZi i2 CZ 2i 3|c + t)i 23) 



(73) 



where the states |c) and |t) could be thought of as control 
and target states respectively for some entangling gate. 
Now, qub its 1 and 3 are measured in the usual Strategy I 
basis (71 1 with £W = 0, while qubit 2 is measured in the 



eigenbasis of the Pauli operator Y, suitably rotated by 
U 1 - 2 ' . This procedure teleports the state initially situated 
on qubits 1 and 3 through an entangling gate 



= R* 

x R r 



(3) 



1 
771 3 



)TSl +7712 



H, 



X™ 2+m3 H 3 M li3 



(74) 



to qubits 4 and 5, with 
M 1>3 = |00> <00| i,3+i |01> <01| 1)3 



-t|10)(10| 



1.3" 



H)(ll|l,3- 

(75) 

Here, the {/jW) are the standard Strategy I byproduct 
angles, Eq. ([34) or (72). This entangling operation is 
related to CZ via 

CZ 1)3 = XrXgR, [7r/2]j R z [tt/2] 3 A/ 1 , 3 X 1 X 3 , (76) 

and so Gij together with single-qubit operators forms a 
universal set of gates. 



C. Probabilistic single-qubit rotations: B — U — B 
state 

Next consider a one-dimensional state of the form 
\R) = B^uf^B^uf • • ^B^+^ICW), (77) 



FIG. 4: B — U — B state, a one-dimensional structure that can 
be used for probabilistic random-length single-qubit rotations. 



where the {B^ 2i+1 ^} are B-type operators, and the 
|U( 2 *)} are once again local unitaries (c.f. Fig. |4|. This 
structure ensures that none of the bonds present in the 
structure is perfect; no particle has maximal entropy of 
entanglement with the rest of the state. This fact is a 
consequence of Lemma |"lII.4| there is no qubit unaffected 
by B-type operators whose neighborhood contains any 
unaffected qubits. 

All single-qubit measurements for the odd-numbered 
qubits now correspond to Strategy II, in which 
the byproduct operator is always X if it oc- 
curs. All even- numbered qubits are measured in the 
{U (2j) |+),U (2l) |-)} basis; as in the previous example, 
the only purpose of these measurements is to enable the 
removal of undesired contributions to the rotation an- 
gles. The main difference from the previous example is 
that a non-unitary gate of the form HR Z [fai+i] is tele- 
ported, where £21+1 G C The present goal is therefore to 
compensate for the imaginary part of the rotation angle. 



As discussed in the previous section and Eq. (60), the 



imaginary part of £21+1 is entirely determined by the ra- 
tio of the singular values of B^ 2l+1 \ and can be defined 
as e = i In (cot #( 2l+1 )) . Consider momentarily the spe- 
cial case where the {B^ 2i+1 ^} all have the same singu- 
lar values. The gate teleported by a measurement of 
the B-transformed qubit 2i + 1 will then be proportional 
to R z [e(— l) m2i +" l2i - 2 + - ], ignoring all rotations about 
real angles which are entirely determined by the choice 
of measurement basis. In short, the sign of the imaginary 
angle depends on the outcomes of the previous measure- 
ments on even-numbered U-transformed qubits. 

The imaginary component therefore undergoes a ran- 
dom walk of step-length |e|. In particular, the walker 
takes its first step to the right when the first measure- 
ment outcome of an even-numbered qubit is 0, and to the 
left if it is 1. Subsequently, a measurement outcome of 
on an even-numbered qubit causes the walker to take 
another step in the same direction as the previous step, 
while outcome 1 makes the walker take a step in the op- 
posite direction. 

The two possible measurement outcomes with odd 
qubits are always equally likely, but the probabilities with 
even qubits depend on the singular values of the B-type 
operators from the (odd) neighboring qubits, and gener- 
ally speaking the walker is more likely to stray further 
from the origin than to step back towards it. For ex- 
ample, consider measuring the first two qubits of \R) in 
Eq. ( [77] ) in the {|+), |— )} basis. It can easily be shown, 
using the Schmidt decomposition ( 38 ) and the expres- 
sion (58), that the probabilities of the outcome |±) on 



qubit 1 are equal, and that those on qubit 2 are 
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p±, 2 = - 1 1 ± cos26» (1) cos26» (3) 



(78) 



Here, the random walk effectively begins at position 
In (tanflW) and moves to In (tan6> (1) ) ±ln (tan6K 3 )). For 
a situation to arise where the walker moves closer to the 
origin with probability greater than 1/2, one of the two 
pairs of conditions 



In (tan6» (1) ) ± In (ti 



> 2 



In 



± cos 26> (1) cos26» (3) > 



(tan# (1 >)|(79) 
(80) 



must be simultaneously satisfied, for either sign. If the 
{0^} are chosen uniformly at random, then the proba- 
bility of this happening is only about 0.315. If this mea- 
surement procedure is continued down the chain, with 
qubits 3-5 relabelled 1-3 after the first two measurements 
and so on, the current value of 6^' tends to drift away 
from 7r/4 towards either or ir/2, and the range of val- 
ues of 9^ for which the walker is likely to turn around 
and walk towards the origin progressively shrinks. Fur- 
thermore, if j^ 1 - 1 } and |#( 3 '} are equal at any time, the 
walker is guaranteed to be more likely to continue in one 
direction than to turn around. 

This procedure constitutes a probabilistic method for 
implementing a single-qubit rotation. Unfortunately if 
the walker strays too far from the origin, it becomes ef- 
fectively impossible to recover and the attempted gate 
teleportation fails. The entire computation must then be 
repeated. If the singular values of the B 2l+1 are chosen 
such that the imaginary components of the teleported 
angles are all integer multiples of each other, then the 
behavior of the random walk is even more deleterious. 
A judicious two-dimensional arrangement of B — U — B 
chains avoids this catastrophe, as discussed in the next 
example. 



D. Universal MBQC: Percolated 2D Cluster State 
from 3D B - U - B state 



In the previous example using Strategy II, a possibly 
infinite number of steps may be required to teleport an 
arbitrary single-qubit unitary. But quitting the proto- 
col results in catastrophic failure: because the computa- 
tional wire is effectively broken, the entire gate teleporta- 
tion must be attempted from the beginning. A solution 
to these problems is to employ the 3D extension to the 
previous resource, corresponding to a cluster state trans- 
formed by alternating B-type operators and unitaries. 
This corresponds to a lattice with two interpenetrating 
cubic sublattices, a B-lattice and a U-lattice. 

An example of this 3D resource, a cube with side length 
3, is depicted in Fig. [5j Initially, Z-basis measurements 
in the z-direction (as labeled in Fig. [5]) are used to carve 
out a structure in which each B-transformed qubit in the 
x — y plane, shaded grey, is attached to a long vertical 





FIG. 5: 3D cluster state with B-type operators and unitaries 
acting on alternate qubits (a) before carving and (b) after 
carving. Grey qubits are acted upon by B-type operators and 
white qubits by local unitaries. Z-basis measurements are 
made in the z direction in (a) to disentangle those vertical 
chains originating from a white qubit in the x — y plane. A 
measurement protocol along the remaining vertical chains in 
(b) produces a percolated 2D cluster in the x — y plane. 



B — U — B chain. Measurements are made on the chain 
qubits, starting at qubit above the B-transformed qubit 
on the computational wire and continuing in the vertical 
direction until success (defined below) is achieved. The 
goal is to probabilistically produce perfect entanglement 
in the x — y plane, thereby effectively eliminating the 
B operators in the horizontal direction. The result is a 
2D cluster state in this plane with missing entanglement 
bonds in random locations. As long as the mean density 
of broken links exceeds the percolation threshold for a 
two-dimensional square lattice, the resource is universal 
for MBQC 

Consider the first vertical chain from the left in Fig. [5] 
Recall that one can interpret B as a ^-rotation by an 
imaginary angle ±zlnA, as shown in Eq. (59). For sim- 



plicity, we assume the unitary operators acting on even- 
numbered qubits are all equal to the identity; were they 
not, they could be compensated by a suitable rotation of 
the measurement basis for. The portion of the state cor- 
responding to the first four qubits of the vertical chain 
(with the qubit that intersects the horizontal chain la- 
beled 1), is then (ignoring normalization) 



\T) = B 1 B 3 CZ 1 , 2 CZ 2i3 CZ 3i4 |+ + + 
= R,pnA] 1 R,[dnA] 3 |Cl 4 )i234, 



'1234 



ignoring normalization factors as usual. First, qubits 2 
and 3 are measured in the {|+), |— )} basis (of course, the 
measurement basis for qubit 2 would need to be rotated 
if a local unitary U^ 2 ) were acting) . These are commuting 
measurements, since they are on different qubits and not 
adaptive. The effect is to teleport the state R z [zlnA] |+) 
from qubit 3 through X m2 HX mi H = X" l2 Z" 11 to qubit 1, 
yielding the new state 

\T') = K z [ilnA]! CZ^xrzrR* [flnAJi \ + 
= CZ 1A R Z [i((-l) m2 + lJlnA)]! | + +)i4 

up to overall local unitaries on the final state. 
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If 7Ti2 = 1, then the imaginary part of the rotation 
angle is completely canceled. Qubit 4 can then be mea- 
sured in the computational basis (again, suitably rotated 
if necessary) to disentangle the rest of the vertical chain 
from the horizontal chain. The result is a perfect cluster 
along the first three qubits in the horizontal direction, 
and the B-type operator is effectively deleted. 

If 7712 = 0, then the situation is similar to the original. 
There is still a B-type operator present in the horizontal 
direction, now corresponding to a z-rotation about an 
angle with imaginary part 21nA = InA 2 . In other words, 
the new effective B-type operator in the horizontal chain 
has a ratio of singular values that is the square of the 
original one. In order to remove the effect of the B-type 
operator, the chain qubits must be measured sequentially 
until the total number of steps towards the origin exceeds 
by 1 the total number of steps away. 

The probabilities and p^ 1 of the outcomes and 1 
on an even qubit in the vertical chain, where k > is the 
present position of the walker on the real number line, 
are given by 



(fe) 
Po = 



1 + \ 2k + 2 



Pi 



(k) 



1 + A 2 + \ 2k + A 2fe + 2 

A 2 + A 2fc 
1 + A 2 + A 2fe + A 2fe + 2 



0(1); (81) 
0(A 2 ). (82) 



The total probability p n that the effect of the B will be 
undone within In measurements is the sum of the prob- 
abilities of all of the possible trajectories of the walker 
on n or fewer steps with initial position 1, final position 
and all intermediate positions strictly positive. 

The probability jjio of undoing the B operator after 
10 attempts (20 measurements) is shown in Fig. [6] as a 
function of the ratio of singular values A. Calculation of 
the exact probability p^ is computationally intractable, 
for two reasons. First, the number of valid trajectories 
for the walker grows exponentially in the number of steps 
allowed. Second, the probability of any particular trajec- 
tory depends on the full history of the walker, not just 
the number of steps. Of course, Poo must approach 1 as 
A approaches unity (the limit that B becomes a unitary 
matrix). In this case, the walk reduces to the simple 
ID random walk, which is known to sample the origin 
frequently. 

If after some predetermined number of measurements 
along a vertical chain one has not yet succeeded in un- 
doing B, the qubit at the root of the chain (i.e. in the 
computational wire) can be measured in the computa- 
tional basis and thereby deleted. The result is a broken 
link in the 2D cluster state. The important result shown 
in Fig. [6] is that there is a critical value of A, called A c , 
above which the probability of successfully undoing the 
B rises above the (bond) percolation threshold for a 2D 
square lattice (approximately 0.593). For this walk, the 
critical value obeys A c < 0.379. The upper bound for A c 
is read off the thick blue curve from Fig. [6j Thus, this 
procedure probabilistically yields a universal resource for 




FIG. 6: (Color online) Probability of success of the procedure 
for deleting a B-type operator in the plane via a random walk 
in the third dimension, as a function of the ratio A of the 
singular values (thick blue, color online). The probabilities 
Pk of deleting B with exactly k even-qubit measurements are 
also shown for k from 1 to 10 (thin, decreasing with increasing 
k), and the thick blue line is the sum of these. The red dashed 
line is the percolation threshold. 



MBQC provided that A is sufficiently large. We note 
that a similar example was considered in |15) . where the 
resource was a 2D cluster state with identical B-type op- 
erators acting everywhere and the percolation proceeded 
via two-element POVMs that either removed the B or 
deleted the qubit. There, the critical value of A was found 
to be 0.649. 



V. DISCUSSION AND CONCLUSIONS 

Motivated by a desire to identify new resource states 
for measurement-based quantum computing, we have 
performed a (non-exhaustive) search of the equivalence 
class of n-qubit cluster states on a rectangular lattice, un- 
der the action of GL (2, C)® n . In particular, our aim was 
to identify which states within this class could be used as 
resources for MBQC, by designing explicit protocols for 
teleporting single-qubit gates and two-qubit entangling 
gates, driven by adaptive local projective measurements. 
We identified a class of one-dimensional states, the so- 
called N — U — N states, that are deterministically uni- 
versal for single-qubit rotations, although with a random 
number of measurements needed to teleport the desired 
rotation. We also identified a probabilistically universal 
resource for single-qubit rotations: the so-called B— U— B 
states. We then described a three-dimensional extension 
of B — U — B states that can yield a universal resource 
for deterministic MBQC beyond a percolation threshold, 
and a 2D N — U — N state that is also universal for de- 
terministic but random-length MBQC. 

Several interesting open issues arise as the result of this 
work. First, it is not clear what (if any) relationship ex- 
ists between the states uncovered in this work and other 
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known resource states. For example, the probabilistic 
nature of the protocol with B — U — B states has fea- 
tures in common with that of other resources for MBQC, 
such as photonic cluster states prepared via probabilis- 
tic entangling gates or with unreliable sources [33H36] . 
Likewise, the quasi-deterministic N — U — N states share 
various characteristics with the AKLT-inspired resources 
of Refs. [5J HOI [HI HD] , m particular the exponential spin 
correlations and the repeat-until-success measurement- 
based strategies. One distinction is that only one alter- 
nating sublattice of the N — U — N states exhibits non- 
zero correlation functions, whereas in the AKLT chain, 
every qubit is correlated with every other. Presumably a 
true identification of an AKLT-type resource with a N- 
transformed cluster state will require N-type operators 
to be present on every qubit, a case we have not handled 
here. 

Also, it will be important to better understand the 
relationships of these states with the universal quantum 
wires of Ref. [T5] . In that work, certain reasonable physi- 
cal assumptions were imposed on ID wires at the outset, 
for example: the possibility of producing a wire via a 
translationally invariant nearest-neighbour global entan- 
gling operation, the asymptotic sharing of an ebit of en- 
tanglement between the left and right halves of the chain, 
etc. In our work, it is not clear if a translationally invari- 
ant scheme exists for producing N — U — NorB — U — B 
chains. Evidence from exact calculations on small chains 
and the explicit description of the states within the MPS 
representation reveals that although the left and right 
halves of N — U — N chains share an ebit, the halves of 
the B — U — B chains do not. This seems reasonable, as 
the N — U — N chain is quasideterministically universal 
for single-qubit rotations, while the B — U — B chain is 
only probabilistically so. 

Second, it is conceivable that all states that have been 
hitherto identified as universal resources for MBQC are 
in fact SLOCC-equivalent to the family of cluster states. 
There is some evidence to support this conjecture. For 
example, the results of Ref. [19] show that many seem- 
ingly diverse resource states can be reduced to cluster 
states via local strategies. Similarly, the proof of the uni- 
versality of the 2D AKLT state on a honeycomb lattice 
proceeds via local reduction to a random graph state, 
which can in turn be reduced to a percolated cluster 
state [13] . This reduction is successful despite the fact 



that the initial resource is defined on qutrits rather than 
qubits, and on a non-rectangular lattice. An even more 
intriguing possibility (though we believe it to be unlikely) 
is that all possible states for universal MBQC fall within 
the SLOCC-equivalence class of the cluster states. At 
the very least, the relative size of this class decreases ex- 
ponentially with the total number of physical qubits, as 
expected [TBI ITT] . 

Third, while we have shown that a certain sub- 
set of the orbit of the cluster states under SLOCC 
are useful resources, either probabilistically or quasi- 
deterministically, it is not clear if the remaining SLOCC- 
transformed cluster states are also universal resources 
for MBQC. Generically, single-qubit measurements on 
these states teleport gates with byproduct operators that 
are rotations about the X axis by complex angles. One 
possibility is that there is a measurement protocol that 
can accommodate all possible byproducts that we sim- 
ply haven't found. Perhaps there is another sense, be- 
sides quasi-deterministic or probabilistic, in which these 
states can be said to be useful for MBQC. Another possi- 
bility is that there is some map from the full orbit to the 
particular subset of states considered in the work. Al- 
ternatively, the states in the full orbit may not be useful 
resources, though we have not attempted to prove this. 

Finally, it would be useful if the resources presented 
in this work could be realized as the ground states of a 
physical Hamiltonian. A recent no-go theorem [37 shows 
this cannot be the case for frustration-free Hamiltoni- 
ans on qubits. The question remains open for frustrated 
Hamiltonians, for instance the AKLT Hamiltonian in the 
Haldane phase [2Tfl - [2"2"] . Another possibility is that the 
ground states of physical Hamiltonians could be locally 
reduced to the resources we have found. Further research 
is needed to answer these and related questions. 
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